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Abstract. We prove asymptotic normality of the distributions denned 
by the central string functions and the basic specialization of fusion 
modules of the current algebra of sb- The limit is taken over linearly 
scaled fusion powers of a fixed collection of irreducible representations. 
This includes as special instances all Demazure modules of the associated 
affine Kac-Moody algebra. Our results are derived from the interpretation 
of fermionic expressions for the characters of fusion products as mixtures 
of probability generating functions. In particular, they describe the 
position of dominating weights and are accompanied by conjectured 
local central limit theorems giving the first hint towards the asymptotic 
growth of dominating weight multiplicities in fusion modules found in 
the literature. Along with an available complementary result on the 
asymptotic normality of the basic specialization of graded tensors of the 
type A standard representation, our result is a central limit theorem for 
a serious class of graded tensors. It therefore serves as an indication 
towards universal behavior: The central string functions and the basic 
specialization of fusion and, in particular, Demazure modules behave 
asymptotically normal, as the number of fusions scale linearly in an 
asymptotic parameter, N say. We remark on the geometric interpretation 
in terms of rational points of unipotent partial flag manifolds over finite 
fields at the end. 



1. Introduction 

The g-supernomials coefficients are ^-analogues of the coefficients of the 
monomials x a in the expansion 

of n^liC 1 + x + ... + x j ) L K They have 
a combinatorial interpretation in terms of generalized Durfee dissection 
partitions in that they are generating functions of so-called (Li, . . . , in- 
admissible partitions with exactly a parts |35| §3]. In there most simple 
form, when m = 1, they are the well-known g-binomial coefficients which 
are the generating functions for restricted partitions [2, §3], a result that 
goes back to Gauss and coins them Gaussian binomial coefficients. Let us 
mention that, geometrically, the coefficients of g-binomials count the number 
of invariants of binary forms [39] §3.3], [U Corollary]. 
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The main focus on g-supernomials, from our point of view, lies on their 
appearance as characters of finite-dimensional modules in the representation 
theory of infinite-dimensional Lie algebras. Precisely, they describe the string 
functions in fusion modules, i.e. graded tensor products, of the current algebra 
SI2 <8> C[t] [121 E3]. Consequently, their coefficients encode the dimensions 
of weight spaces, i.e. isotropic components with respect to the action of a 
maximal abelian subalgebra. 

The exact determination of those coefficients is certainly possible in any 
fixed instance of a g-supernomial, simply by extrapolation from the definition 
of g-supernomials and their building blocks, the g-binomials. In general 
though, the explicit description of those coefficients remains intractable and 
one usually is satisfied with concrete expressions for their generating function, 
the g-supernomial. We will examine those coefficients from a qualitative point 
of view by the interpretation of g-supernomials as probability generating 
functions of discrete distributions. Furthermore, we investigate the gener- 
ating function of the q-supernomials themselves. This generating function 
equals the so-called basic specialization of the character of SI2 <8> C[t]-fusion 
modules, and can be understood as the Hilbert polynomial associated to the 



parts that are graded by the action of t. Our main results are Theorem 4.20 



and Theorem 4.24, and can be summarized as follows: 



Consider a sequence (C 1 )* 6 ^ * (C 2 )* L i JV) * • • • * (c m+1 )* L ^° of fusion 
modules of the current algebra 5I2 <S> C[t]. Assume that the exponents grow on 

a linear scale, i.e. {l\ , . . . ,L^)/N ->a/0. Then, the sequences associ- 
ated to the central string functions and basic specialization of those modules 
behave asymptotically normal with mean and variance growing quadratically 
and cubically in N, respectively, with explicitly calculable leading terms. 

Along with complementary results for the Demazure module V-n Ui (Aq) of 
the higher rank affine Kac- Moody algebra si r [3 §5.4], our findings highlight 
the first serious indication of a central limit theorem concerning the central 
string functions and the basic specialization of fusion modules. We will 
conclude this line of thought by a discussion of fusion of symmetric power 
representations in type A and the following conjecture. 



Conjecture 6.5 , Consider the sequence of type A fusion modules of sym- 



metric power representations 



„AN) f (N) AN) 



IN) f W AN) AN). (j\n 

associated to the partition [P- > = (l/i ,2^ ,...,&h ) with f,- -many 
i 's. Assume that as N — > 00 we have 
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Then, the central string functions and the basic specialization of F^n) behaves 
asymptotically normal as N — » oo. 

Since Demazure modules of sh are special instances of fusion modules we 



will compare (see §4.5.2) our results to previously derived formulas for the 
expectation value of the basic specialization of Demazure modules |7j. There 
the expectation value has been derived via different methods, through the 
detailed analysis of the recursion given by Demazure's character formula. 

Our results are accompanied by conjectured local central limit theorems 
in §5, that give the first hint towards the asymptotic growth of dominating 
weight multiplicities in fusion modules found in the literature. We furthermore 
discuss the geometric interpretation in terms of rational points of unipotent 



partial flag manifolds over finite fields in Remark 6.4 



2. Notation 



2.1. g-supernomial coefficients. We follow [35] and fix a vector L = 
(L\, . . . , L m ) G Z™, and let l m = YliLi^i- Consider the expansion [35 
(2.5)] 



Y[(l + x + ... + x j ) 

3=1 



and the fermionic representation |35l (2.8)] of the coefficient 



E 



r 

3m— 1 



Ll + 32 

31 



JiH hj m =a+- 



as a convex sum of products of usual binomial coefficients. Such an expression 
as a positive sum of products of binomials is commonly referred to as fermionic. 
Schilling and Warnaar define the g-supernomials [351 (2-9)] as a g-analogue 
of those coefficients 



(2.1) 



Here j m+1 = 0, 



E 

JlH hjm=a+-f 

Ml 



+ ji+i 
3l 



[M\ q \ 
[k] q \[M-k] q \ 



denotes the q-binomial coefficient and 



[k]q\ = rij=i T=^" the g-factorial (see e.g. [10]). Note that Feigin and Feigin 
describe in |12t Theorem 5.1] the characters of sl2 ® C[t] fusion modules in 
terms of a slight modification of [ Q ] , that is [33 (3.1)] 



h a 



T(L,a)(q)=q L i 
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with its explicit form [35,, (3.2)] 
(2.2) f(L,a)(q) = 



n 



Lt + je+i 
k 



E < 

JlH hjm=a 

Again, we let j m +i = 0. We will discuss characters of fusion modules as 
generating functions of g-supernomials in greater detail later in §3. This 



explicit expression of T can be derived easily from the above definition of 
[^] due to the symmetry of the g-binomial coefficients 



For illustration purposes, in the case m 
give 

\(Li) 
a 



E 



-k) 


~N~ 








k_ 










g 




- 1 the definitions 


~L{ 






Li 


.31. 


q 







f{{L x ),a){q) 



ji=a+- 



~L{ 




~L{ 


Ji. 


q 


a 



That is, T is a shifted (here by L\/2) and translated (here by a 2 ) version of 

r(^i)i . 

L a J q 

2.2. Lie algebras, Demazure and fusion modules. For general facts 
about current and affine Kac-Moody algebras and their representation theory 
we refer the reader to [23]. We denote by s[ r the complex- valued r x r 
matrices with trace (the reader might adapt to the case r = 2 as this 
will be the setting we will mostly consider). Then, sir ® C[t] denotes its 
current algebra and, closely related, sl r its associated affine Kac-Moody 
algebra which can be realized as the extended loop algebra of sl r , i.e. s\ r = 
Sir ® C[M -1 ] ©Cd©Cc (see [231 §7]). Here, d denotes the derivation tf t 
and c the canonical central element. In view of that realization let b D f) be 
the Borel and Cartan subalgebra in s\ r corresponding to their a priori fixed 
counterparts in sl r . We denote the simple roots by cti, • • • , a r -i £ f)*, the 
highest root by 6 = ot\ + . . . + a r -l> the imaginary root by 5 = ao + 9, and 
the simple coroots by otQ , . . . , G f). Note that we assume ai, . . . , a r -i 
to correspond to the standard embedding of s\ r . Let sq, si, . . . , s r -i be 
the simple reflections associated to the simple roots and let the subgroup 
W aS = (s , si, . . . , s r -i) of GL(h*) denote the affine Weyl group. 

For a dominant integral weight A = miAo + r^Ai + . . . + m r _i A r _i in the 
affine weight lattice P = ZA ZAi ... ZA r _i C f)* we let F(A) be the 
integrable highest weight representation of weight A of si r and x(V{A)) its 
character. All weights /x occuring in V(A), that is, having non-trivial weight 
space V(A)u, are elements of the lattice T = A+Zao + Zai-I-- • - + Za r _i C P. 
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One usually writes the monomials in the characters of such modules as e M , 
the coefficient k in the monomial e~ ka ° is referred to as the degree. The 
Ao, Ai, . . . , A r _i are called fundamental weights, the V(Aj) the fundamental 
representations and V(Aq) the basic representation. 

We will not be directly concerned with the representation theory of s[ r 
but instead with Demazure and fusion modules, i.e. finite-dimensional rep- 
resentations of the Borel subalgebra b and the current algebra si r ® C[t], 
respectively (see e.g. E2J Q31 O US 125). 

For w G W aS and a dominant integral weight A = miAo + TO2A1 + 
. . . + m r _iA r _i denote the associated Demazure module by V W (A). Such 
a Demazure module is defined as the finite-dimensional subspace of the 
integrable highest weight representation V(A) generated by the action of 
the universal enveloping algebra U{b) on a non-zero weight vector v w ^ 
inside the 1-dimensional weight space V(A) W ^, in short V W (A) = U{b).v w ^y 
Demazure's character formula [281 131 j allows the computation of the 
character x(V w (A)) by an iterated application of certain operators on the 
monomial e A = x(^i(A)), as follows. We introduce the convention that 
Yl7=o a « = an d Yli=o a i = ~ a -i — • • • — a/c+i for k < —1. Note that this is 
natural in the sense that Gauss's summation formula ^^=0 * = extends 
to all k G Z, as does the identity Y2i=o 1 = k + 1. With this convention 
the Demazure operator Dj associated with a simple reflection Sj acts on 
monomials as 

(2.3) D je ^ = e>*- ia i. 

i=0 

For an arbitrary Weyl group element w G VF aff we choose a reduced decompo- 
sition w = Sj 1 Sj 2 ■ ■ ■ Sj t and set D w = Dj 1 Dj 2 • • • Dj r Demazure's character 
formula now states that the character of V W (A) can be computed recursively 
as X (V w (A)) = D w e A . 

For a general definition of a fusion module consider irreducible st r -modules 
V(Ai), . . . , V(\ n ) with integral dominant weights A = (Ai, . . . , A n ). For 
regular Z = (zi, • • • , z n ) G C n let sl r ® C[t] act by 

(x <g) t r ).v = zl(x.v), for x G si r , v G V(Aj). 

This so-called evaluation module is denoted by ^(A^). Then \\{Z) = 
n) is a cyclic module of the universal enveloping algebra 
U(slr ® C[i\). Hence, there is an induced filtration 

V A (Z) s = 0^(sl r ®C[%v 

3=0 

The associated graded module is called the fusion module: 

F X {Z) = t^(Ai) * • • • * ^ n (A n ) = grV A (Z) = V X (Z) S /V X (Z) S . 1 . 

s>0 



(i 
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The additional grading is induced by the powers of t. We obtain a graded 
tensor product multiplicity £ Z[q] such that 

where m M is the usual tensor product multiplicity of V([i) in <S>iY(Aj). For 
specific A = (Ai, . . . , A n ) it is known that the fusion modules J-\(Z) are 
independent of the choice of the evaluation parameter Z, see e.g. §1.5.1], 

2.3. Statistical notions. Standard sources are [3[Hl ll5U16| . All our random 
variables X will be discrete and finite. Recall that the expected value of 
such a random variable is the weighted average E(X) = ^^{X = x)x. 
The covariance of two random variables X and Y is Cov(X, Y) = E((X — 
E(X))(Y - E(Y))). They are said to be uncorrelated if Cov(X, Y) = 0. The 
variance of X is Var(X) = Cov(X, X). Its probability generating function 
is E(q x ) = ^ X P(-X" = x)q x , and the associated probability distribution 
fix = Ylx^i-^ = x )$x- Here, 6 X denotes the Dirac distribution (point mass) 
at x. A sequence Xx converges P-almost surely (P — a.s. for short) to X 
if P(lirn/v-).oo Xx = X) = 1. Convergence and equality in distribution will 

be denoted by — > and =, respectively. Af(fi, E) will denote the normal 
distribution with mean n and covariance matrix S. Note that Af(fi, 0) = 5^. 
The conditional probability P(Y = y\X = x) = P(X = x)- l ~P(X = x, Y = 
y) is the probability of Y taking the value y given the occurence of the value x 
for X. A mixture of distributions \iXi is a finite convex combination thereof, 
i.e. ^ZiWiUXi for some weights W{ > with Yli w i = 1- The probability 
generating function of a mixture is Yli wiE{q *)■ 



3. Characters and their specializations 

3.1. Graded characters of fusion and Demazure modules of s^. We 
focus on [12l [13] which study in detail fusion modules for the current algebra 
eh <8) C[t]. Note that irreducible representations of are indexed by natural 
numbers No = {0, 1, 2, 3, . . .}. Hence, a fusion module F&(Z) for 5(2 is given 
by a dimension vector d = (di, ■ ■ ■ , d n ) 6 Nq and regular Z £ C n . For sl% it 
is known that the fusion modules are independent of the choice of Z € C n 
|12j . and it is usual to omit the evaluation parameter in the notation. In 
particular, Feigin and Feigin [12] study the fusion module J-^ with dimension 
vector d = (I 61 , 2 te , . . . , (m + l) 6 "^ 1 ), and denote it 

T(b 1} 63, ... , 6 m+1 ) = (C 1 )*' 1 * (C 2 )*' 2 * • • • * ( C m+1 )* Wl . 
They prove |12l Theorem 5.1] that, written here in terms of the q-supernomial 



T(L, a)(q) (2.2), its graded character can be described as 



(3.1) x(Hh,b 2 , . . . , b m+l ))(z, 5 ) = Vz". T(L, a)(g). 

V- s £ * 
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It is well-known that Demazure modules V W (A) associated to SI2 carry a 
sl2<£>C[i]-module structure and as such are special instances of fusion modules 
(see e.g. [HI §1.5.1] or [HI §3.5]). To be precise, we have isomorphisms of 
sfe <K> C[t] -modules as follows 



(3.2) K,(mA + nAi) 



J"(0, 1 ,0,l(w) - 1) ,w = w's , 

m + 1 m+n+l 

J"(0, 1 ,0,l(w) - 1) ,w = w'sx. 

n+l m+n+l 



Here, we write l(w) for the length of a reduced decomposition of w G W . 
All elements of W aS have the form 

wn,q = • • • sosiso or w N .i = • • • SiS Si 

V v ' V v ' 

N factors N factors 

for N >Q. Hence, Demazure's operator D w becomes an alternating iteration 
...DqD\DqD\... of the operators Dq,D\ defined in (2.3). Let us summarize 
the non-trivial elements of the affine Weyl group W aS of sh as 

( Sl s ) N Sl ,N > and (s oSl ) N , (s lSQ ) N , s (s 1 s ) N - 1 , N > 0. 

Then, one can precisely identify the characters of Demazure modules via 
those of fusion modules as follows. 

Proposition 3.1. Consider the Demazure module V w = V w (mAo + nAi) 
of SI2 of (fixed) highest weight A = itiAq + nA\ and recall the character 
formula for fusion modules (3.1). The character x(Vw)(z,q), written in the 
coordinates z = e~ ai , q = e~ d , is given by 

(3.3) For the trivial element w = 1 one has x{Vi){ z iQ) = 

(3.4) For w = (siso) si, N > one has 

e- A X (V w )(z,q) = 

2 -(n+ m) JV-n/2 g ^ m+ iV(i V+1 ) Tlx(;F(0) Lw ^ zq -2N-l^ q) 

where L w = (Li, . . . , L m+n ) = (0, 1,0, 2N). 

n m+n 

(3.5) For w = (sqs{) n , N > one has 

e' A X (V w )(z,q) = 

z - { n+m)N-n/2 q N*m+N(N+l)n x (jr^ L ^ ( zq -2N-l ^ ?) 

where L w = (Li, . . .,L m+n ) = (0, 1,0, 2N - 1). 

n m+n 

(3.6) For w = {s\Sq) n , N > one has 

e- A x(V w )(z,q) = 

z -(n+m)N+n/2 q N*m+N(N-l)n x (jr^ Lt|) )( zg -2W q) 

where L w = (Li, . . . ,L m+n ) = (0, 1,0, 2N - 1). 

m m+n 

(3.7) For w = sq(siSq) n ~ 1 , N > one has 



s 
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- A x(V w )(z,q) = 

z -(n+m)N+n/2 q N*m + N(N-l)n x (j:^ Lw )( zq -2N ^ g) 



where L, 



(Li, . . . , L r 



(0, l,0,2iV-2) 

m m+n 



The sum of the entries in L w represents the length l(w) of the Weyl group 
element w. When either n or m equals 0, then L w = (0, l(w)). 

Proof. Feigin [13, (11)] denotes an integrable highest weight representation 
Li t k = U(sl2)-Vi ! k with highest weight vector v^k such that c.v^k = kvi t k, 
ho- v i,k = ivi^k, and d.Vi^ = 0. In our notation, the canonical central element 
is c = Oq + a\ , the coroot is ho = Oq , and the scaling element d is given by 
ao(d) = 1 and ai(d) = 0. Therefore, by comparison of the highest weight 
vector we have L,^ = V((k — i)Ao + iAi). The bigrading is chosen according 
to the action of ho and d, and consequently, the character is denoted in the 
monomials e a ° = e s ~ ai and e _<5 , respectively. By |13l Corollary 3.1] each 
such module can be constructed as an inductive limit of fusion products, 
that is Lj^ = C t+l * (C fc+1 ) 2o °. Each fusion product can be identified with 
the corresponding Demazure module V w ((k — i)Ao + iAi) by comparing the 
weights of the extremal weight vectors described in |13|, § 1] . Now apply the 
character formula \12\ Theorem 5.1], noting that e a ° = zq~ l . □ 

Let us illustrate g-supernomials, characters of Demazure and fusion mod- 
ules, and their relations by some examples. We consider m = 2, n = and 
the vector L = {L\,L2) = (0,4). Let us compare: 



^1) The generating function of the g-supernomials [^ 4 ^ 



f L (z,q) = Y,z a - 



aeZ 



(0,4) 
a 



Note the explicit form (compare (2.1)) 
~ (M) 1 = V <t 

q Ji+j2=a+4 



32) 


"4" 




~32 




h_ 


q 





and see Figure 1 for the actual plot. 



(2) The character of the fusion module ^-"(0, 0, 4) = (C 3 )* 4 as depicted 
in (3.1 ), i.e. 

X(.F(0, 0, 4))(z, q) = £ z a ■ T((0, 4), a)(q). 

aeZ 



Note the explicit form (compare (2.2)) 



r((0,4),a)(g)= Y, 



"4" 




~32 


J2_ 


9 


Ji. 



and see Figure 2 for the actual plot. 
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Figure 1. Plot of the generating function of q-supernomials fi,(z, i 



(3) The character of the Demazure module V( SlSo \2(2Ao) written as a 
Laurent polynomial in the monomials z = e~ ai ,q = e~ s , i.e. the 
generating function of its weight space dimensions, 

x(V {siSo) 2(2A ))(z,q) = *V • dim (V (siSo)2 (2A ) z a q i) . 

a,i£Z 

This generating function can be compute via Demazure's ( recursive) 
character formula x(V( sl s ) 2 (^o)) = DiDoD\Doe 2Ao . See Figure 3 
for the actual plot. 
As to the relations between the different generating functions, one has 

(3.8) e- A X (V (siSo) 2(2A ))(z,q) = q 8 • f L (z,l/q), 

(3.9) e- A x(V (siSo)2 (2A ))(z,(z) = z-V-x(^(0,0,4))(zg- 4 ,(?). 

All expressions are equivalent up to translations (multiplication by e.g. z~ 4 ), 
reflections (evaluation at the reciprocal 1/q), and rotations (evaluation at 
mixed monomials zq~ 4 ), respectively. Note that the example (3.9) is an 



instance of the general formula (3.6). 



3.2. Real characters and the basic specialization. As already men- 
tioned in §2.2 the specialization of the character x{F(b\, . . . , b m )) of a fusion 
module F(b u ...,b m ) = (C 1 )^ 1 * (C 2 )* 62 * • • • * (C m+1 )* fe ™+i a t g = 1 equals 
the character of the tensor product (C 1 )® 61 <g> (C 2 )® 62 <g> • • • <g> (C m+l )® b ™+ 1 
of irreducible representations of sfc. Note that frequently this type of spe- 
cialization is called a real character 



(3.10) 



Therefore, (3.2) implies 

x(V w (mA + nAi))(z,l) 

Jx(C m+1 (8) (C m+n+1 



x(c 



n+l 



/^m+n+lj®I(tc)-lw^ 



,W 

, w 



w s , 
w'si. 
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Figure 2. Plot of the character x(.F(0, 0, 4)). 

This so-called factorization phenomenon for Demazure modules is known 
independent of the theory of fusion modules, first proved in [32], 33J, and 
generalized considerably, e.g. |18j . Note, since characters of (non-graded) 
tensor products multiply, i.e. x(Y ® W) = x(Y) ' x(W), the associated 
distributions convolute and asymptotic considerations, as l(w) — > oo, take 
the simplest form of a central limit theorem for sums of i.i.d. random variables. 
Therefore, characters of usual tensor products, and equivalently the above 
specializations at q = 1, are well understood from a statistical point of view 
and have been analyzed further in great detail, e.g. |42| . 

Much less studied is the so-called basic specialization of those characters, 
i.e. their evaluation at z = 1. We borrow this terminology from Kac |23^ 
§1.5, 10.8, 12.2] who analyzed this kind of specialization for characters of 
integrable highest weight modules ^(A), and obtained Macdonald's identities 
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FIGURE 3. Plot of the character x(V( Sl s ) 2 ( 2j V))) written in 
the coordinates z = e~ ai and q = e~ s . 



for Dedekind's ^-function [23\ §12.2]. Note that the evaluation of Demazure 
characters at z = 1 can be equivalently described as their restriction to the 
subspace Cd when viewed as functions on b = b © sl r © iC[i] © Cd © Cc, 
in short K„(A)| C(J S N[g = e _<5 ]. For specific Demazure modules the basic 
specialization has a geometric interpretation in terms of Galois numbers. That 
is, consider the translations t-jv Wl = (S1S2 • • • s r _i<7 r-1 ) Ar in the extended 
affine Weyl group of sir, where u)% = Ai — Ao and a denotes the automorphism 
of the Dynkin diagram of sl r which sends to 1. We denote the Demazure 
module associated to the fundamental weight Ao and those translations by 
^-iVwi(Ao)- Then, the basic specialization x{V-Nw\ (Ao))|cd ^ N[q = e~ s ] 
equals the so-called generalized Galois number that counts the number of 
flags = Vq C V\ C • • • C V r = of length r in an N- dimensional vector 
space over a field with q elements [27], 143] . For r = 2, that is in the 
case of Demazure modules of sfe, this is the Galois number that counts the 
number of subspaces in finite vector spaces that has been studied initially 
by Goldman and Rota [20]. In general, due to a result of Shimomura [37] . 
(/-supernomials count the number of points of unipotent partial flag varieties 
defined over finite fields (see §6 and, in particular, Remark 6. 4| therein). 



4. Distributions defined by (7-supernomials 
Let L := (Li, . . . , L m ) € and consider the (/-supernormal given in 



(2.2): 



f (L, o)(g) = Ll) II 

jlH Yjm=a. £=1 



Le + Ji+i 
H 
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Recall that these generating functions were introduced by Schilling and 
Warnaar [35J who showed that they enumerate L-admissible partitions with 
exactly a parts. In this section we study the average behavior of the distri- 
butions defined by T(L,a)(q) and by the generating function 

Li+...+L m 

f(L)(q):= Yl nL,a)(q) 
of the total number of L-admissible partitions that equals the specialization 



x(-7~~(0, L))(l, q) by (3.1). We show that the total number and in certain 
(typical) cases the a-restricted number of L-admissible partitions are asymp- 
totically normally distributed with asymptotic parameter being a convergent 



sequence ^L^, as N — > oo. 

4.1. Preliminaries. The distributions with probability generating function 



F a ,b{q) ■= [ a a b ] l{ a t b ) were first investigated by Mann and Whitney [30], 



who showed: 

Theorem 4.1. Let Inv a ^ be a random variable with distribution F a ^. Then 
Inv a j, has expectation ¥j{Inv aj b) = ^ab, variance Yai(Inv a fi) = ^a6(a + 6 + 
1), and 

Inv n h — ~E(Inv n b) d 
y/Vax(lnv a , b ) 

as a, b —7- oo. 

A corresponding local limit theorem was proved by Takacs [41] . 

Remark 4.2. The (/-binomials enumerate different number theoretic [2], 
geometric [39], and combinatorial [JT] objects. We interpret them here 
as counting inversions |4U) . Consider a word (unordered sequence) w = 
(wi, . . . ,w n ) of elements from an ordered set. A 4-tuple (i, j,Wi,Wj) with 
i < j and Wi > Wj is called inversion. It is well known that [ a ^ fe ] 9 is the 
generating function for inversions in words of a zeroes and b ones. 

We call a vector j = (ji, . . . , j m ) £ NJJ 1 compatible to L G if ji < 
Li + jj+i for i = l,...m, and a probability distribution L-compatible if 
the set of L-compatible values has probability one. For L-compatible j we 
let Im>(L, j) denote a random variable with probability generating function 
F(L,j)(q) := J]^Li F Li+j i+u ji(q) (here and in the sequel j m+1 = 0), and we 
let 

m i—1 

Q(L,j) = Y, ji { ji + Y, Le )- 

i=l 1=1 

We view the normalized fermionic expression g{q) = T(L)(g)/T(L)(l), and 
g a (q) = T(L, a)(q) /T(L, a)(l) , respectively, as mixtures of probability gener- 
ating functions, i.e. 

5(«) = E^ (LJ)i? ( L 'j)^) p ( J= j)' 
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weighted by 

P(J = j) = P( Jt =j 1 ,...,J m =j m ) = f[ ( Li + Ji+1 ) /f (L)(l). 

All generating functions and random variables considered here depend im- 
plicitly on the admission vector L, but this dependence will from now on for 
convenience be suppressed in the notation. We may equivalently write g{q) 
as follows. 

Proposition 4.3. With given conditional distribution P(Y = i | J = j) = 
P(Inv(L, j) = i) we have 

(4.1) g(q) = B(q^ +Y ). 

Proof. By the definition of the conditional distribution 

5(g) = E^ (L,j)F ( L >j)^) p ( J= j) 

j 

= ^( L J)E(/|J=j)P(J=j) 
j 

= q Q ^ +y P(F = y | J = j) P(J = j) 

j,y 

= Y^ q Q^+yp(Y = y ,j=j) 

j,y 

= E(g^ L ' J ) +y ) □ 

Our interest lies in the distribution of the random variable 

T = Q(L,J)+Y. 

We write Y = E(y|J) + R, where R = Y — E(Y| J) is (by the properties of 
conditional expection) uncorrelated to (any square-integrable function of) J. 
Furthermore, by Theorem 4.1 we have 

e(L, J) = E(Y\3) = X - (^(L m - J m )J m + ( L i + Ji+i ~ Ji)Ji \ ■ 

Let us call the distribution of J the mixing distribution (for the total number 
of L-admissible partitions) . To understand the behavior of mixing distribu- 
tions and to introduce a convenient way to refer to them, let us describe a 
simple random experiment. We focus our attention on the unrestricted case 
first. 

4.1.1. The basic probabilistic setup. Consider m mutually independent ran- 
dom sources Si,S%, . . . , S m emitting words W(l), . . . , W(m). Each word 
W(i) = (Xi(i), X2(i), ■ ■ ■ ,XL { (i)) is a sequence of L« many mutually inde- 
pendent letters Xk(i), where each X^(i) is uniformly distributed from the 
alphabet {0, . . . , i}. 
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For < i < k let B{(k) := ^2j=i denote the random vari- 

able the number of appearances of letter i in word W(k). Then, the ran- 
dom vector B(fc) := (Bo(k), . . . , B/,(k)) is the occupancy (statistic) of word 
W(k), S{k) := J2i=i x i( k ) = T,j=aj B j(k) is the sum of word W(k), 
Bl := (B(l), . . . , B(m)) is the total occupancy of W(l), . . . , W(m), and 
5l := X^=i Sty) 1S called the total sum of words W(l), . . . , W(m). Clearly, 
under the assumptions above, <5x is the sum of independent uniformly dis- 
tributed random variables and we have 



1 m 

(4-2) E(S L ) = -^iLi, 

i=l 

j m 

(4.3) Var(5 L ) = -^(i + 2)^. 



12 

i=l 



Since fi(x) = (1 + X + . . .+x l )/ty + l) is the probability generating function 
for the uniform distribution on {0, . . . , i}, the coefficient of x a in the product 
Y\i=i(fi( x )) Li gives the probability that 5l = a, which we write as 



m 

\.r"] 

i=l 



a ]Y[(fi(x)) Li =P(Sj, = a) 



Now, Schilling and Warnaar [65\ (1.3)] define supernomial coefficients ( L ) 
(for a + i m /2 6 Z+) by 



m oo \ 

n(i+*+...+*o £ '= e C) 

i=l a=-i m /2 V 7 



x a+i m /2 



where ^ m = Yli=i i^i = 2E(5l). We may thus interpret the supernomial 
coefficients as the probabilities of the random variable <Sl : 



)/(n( i + 1 ) Ll ) = P(^L-E(5 L )=a) 



Schilling and Warnaar [351 (2-8)] give the representation 

T \ m /T i • 



+Jm=a «=1 



which leads to 



\ii+-+im=a 4=1 V ^ 7 / \i=l 

and the following result. 



DISTRIBUTIONS DEFINED BY Q-SUPERNOMIALS 



15 



Lemma 4.4. For j = (ji, ■ ■ ■ , j m ) £ N™ (Wid u>i£/i the usual convention 
about binomial coefficients that (f) = unless < b < a) let 

p( Ji = ji, . . . , j m = jm) = (n ( L * +.f +i ) 1 / (rp + • 

Then, these numbers define a L- compatible probability distribution. 
Proof. Indeed, we have a probability distribution since 

Li + ...+L m 

P(Ji=ji,...,J m =j m )= £ P(S L = o) =1. 



The L-compatibility is obvious. □ 
We may describe this probability distribution alternatively as follows. 

Proposition 4.5. Let Bl = (B(l), . . . , B(m)) be as above. That is, the 
random vectors B(i) = (Bq(i), . . . , Bi(i)) are independent, where each B(z) 
/jas a f uniform ) multinomial distribution with parameters Li and po = ■ ■ ■ = 
Pi = jlj. Let 

m i 

(4.4) ^ := 5^ fe _ i+1 (fc), w/iere := J^fij(i)- 



Then, the joint distribution of (J\, . . . , J m ) is as m Lemma 4-4 



Proof. We use (formal) generating functions. It is clear that the joint 
generating function of (Jj, . . . , J m ) as defined in (4.4) is 

m i i m 

m? -e)=(n(E n ^) Li )/(n^+ i ) L,; )- 

i=l j=0 k=i— j+1 i=l 

Now extract coefficients to see that this corresponds to the distribution 
defined in lLemma 4.41 □ 

Remark 4.6. Since A k -i+i(k) = B k ^i + i(k) + B k -i +2 {k)-\ \-B k (k) counts 

the number of appearances of the highest i letters in word W(k), the Ji may 
be described as the total number (overall count) of the i highest non-zero 
letters in all words. 

[Proposition 4.5 shows that the mixing distribution for the total number of 
L-admissible partitions may be realized as a simple linear transformation of 
Bl. We call Bl the underlying occupancy distribution. This representation 
can be used for explicit calculations, and reduces the asymptotic treatment 
of J in the unrestricted case to the well known asymptotics of multinomial 
distributions. Let us recall the following classical result about the asymptotic 
normality of multinomial distributions. 
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Theorem 4.7. Let have the multinomial distribution with parameters 

N and p = (po,Pi, ■ ■ ■ ,Pm)- Then, we have mean E(B^) = Np, covariance 
matrix Cov(B^)) = NT,, and 

BW-ATp d 



ATl/2 

where £ = diag(p) — p*p. 



AA(0,S), 



Given these initial observations, we have a straightforward program to 
treat the asymptotics of the random variable T = Q(L, J) +Y: 

(1) split T into an "occupancy part" Q(L, 3)+E(Y\ J), dependent only on 
J and a remaining "rest-inversion part" = Y — E(Y|J) "orthogonal" 
to J, 



(2) find the asymptotics of the two parts, using Theorem 4.1 and Theo- 



rem 



T7\ 



(3) combine the results. 



This program is carried through in §4.2 



4.1.2. The probabilistic setup for the a-restricted case. It is clear that the 
same experiment describes the a-restricted case when we consider only the 
outcomes with total sum 5l = a. That is, for the a-restricted case the 
underlying occupancy distribution is Bl \ Sj, = a, i.e. the distribution of 
Bl, conditioned to have Sl = a. In order to have a succinct wording for the 
"most important" restricted cases we make the following definition (Cf. (4.2)). 

Definition 4.8. We call the a-restricted cases with a = E(5l) (resp. a = 
E(5 L ) ± \) central. 

This terminology is justified on the one hand because the distribution of Sl 
is symmetric around E(5l), and on the other hand, by the strong law of 
large numbers, we have P — a.s. as h^ N ' — > (ai, • • • , a m ): 



N u 2 

k=l 

So these cases are of central importance. 

4.2. General asymptotic considerations. Throughout this section let 
be a sequence of admission vectors, and a sequence of L^- 

compatible mixing distributions. We consider the sequence of inversion 
statistics 

y (iV) , defined by P(7^ = i | = j) = P(Inv(L^ N \ j) = i), 



and recall the associated definitions from §4.1, in particular ( |4.l| , 

R (N) = yW_ E(F (iv)| J (JV) )) 

TW=Q(LW,jW)+yW, 
e(L W 7 jW ) = E(y W| J W ) . 
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We first show a result that enables us to treat the occurring quadratic 
functions of J^. For quadratic functions of asymptotically normal random 
vectors one has in general: 

Proposition 4.9. Assume there exists b G R m and a positive semidefinite 
matrix S G p^ mxm of positive rank such that 

E(XW) , , XW-ATb ,„ , 

Let M G R mxm , v G R m , anc? consider the quadratic function q(x,v) = 
x*Mx+vx*. Assume additionally thatv^ is a sequence such that -hv( N ' — > 
a and let w = b(M + M*) + a. Then, 

- — 1 — 4^ — - — - A ^°' wSw )' 

Proof. Let 

9 (xW,vW)-g(E(X( ff )),v( w ))=#)+#), 

where 

AW = (E(XW)(M + M*) + vW) (XW _ E (X( 7V )))*, 
BW = (XW - E(#))) M (XW - E(XW))*. 

Now, clearly ^ «, f (jV(0,E))i = AA(0,wSw*) and fg- -4 0, 

and the claim follows. □ 

For the convergence of moments we have here: 

Proposition 4.10. In the situation of \Proposition assume that addi- 
tionally ^Cov(X^) — > S. Then, 

ana 7 , 

E(,(xW,vW)-,(E(xW),vW)) ^ 

If furthermore E (x^ - E(xf°)) /N 3 — > /or all i, then 

Var(g(xW,vW)) t 
„ r WLW . 

We omit the elementary proof, and now look at the asymptotic behavior 
of the conditional distribution of R( N ' given that J^) =j( N ). 

Lemma 4.11. If iL^ — > a ^ 0, andj^ is a sequence of "L^ N ~> -compatible 
vectors such that — >• b 7^ 0, then 

/ra(LW,jW) _ E(I™(LW,jW) d 
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where v(a, b) = ^ YT=i( a i + ~~ h)h{ai + h+i) (where b m+ i = 0). 
Proof. By definition, the random variable 

/™(LW,j( N ))-E(/™(LW,jW)) 
is distributed like the sum YliLi ^1 of m independent random variables 

Let o| = —Ji and = . By assumption the limits dj 

lim ^fl^ and 6j = lim jjb^ exist. If dj = or b-i = 0, then J 5 ^^^-* - 
= AA(0,0). If oi > and 6, > we find that ^Var(/nw(aJ 7V) , b\ N) ) - 



ka>ibi(ai + bi) := w(cii, 6j), and Theorem 4.1 gives that 



12 



N(0,w(ai,bi)). 



N 3 / 2 

Thus under the conditions above 

A/"(0,?j(a,b)). □ 

Remark 4.12. The case Cjdj = for all 2 (that is v(c, d) = 0) is less 
interesting but not excluded. In this case we interpret A/"(0, 0) := 5o as the 
Dirac-measure (point mass) at 0. 

Finally, we turn to the combination of the results above. It turns out that 
under mild conditions the limiting distributions of the (normalized) random 
variables R^ N ' and J^-* are asymptotically independent. 

Theorem 4.13. Let R(&, b) denote a random variable with distribution 
W(0, v(a, b)). // ->a/0, and if there exists b G R + m and a 

positive semidefinite matrix X E f^ mxm of positive rank such that 

^( > S )" 

Then, as N — 00, 

(f5, J ^)^W^(°.s>), 

where the constituents on the right-hand side are independent. 

Proof. Let i C R be a Borel set and / : R m — > R be bounded and 
continuous. We have 



_ E p / fa(LW , JW) - e(LW , JW) g A \ fjW - Nh 



N 3 / 2 J \ N 1 / 2 
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By Skorokhod's representation theorem (see jHEEj) we may assume that 
J( ^ ) 1 7f b — > X, P - a.s., where X = Af(0, S). Then clearly ijW b, 
P — a.s., and by the preceding lemma 

P ( ^ ! ? >- ei)^ A/\0, u(a, b)) (A) , P - a. S . 

Therefore, 

/ R( N ) \ / tW — /Vh\ 
E1 ^y/(^) ^AA(0,,(a,b))(A)E(/(X)). □ 

Let us collect some immediate corollaries. 

Corollary 4.14. In the situation of Theorem 4-lS\ assume that additionally 

^E(JW) b. Let b = b m+1 = 0, St = 0, ^ = |E,=i + 

E(J^) — E^ 7 ^)) E(j£ ), and c = c(a, b) 6e £/ie vector with coordinates 
d = ai + + - 2bi. Then, 

(45) ^.^-^ ^^,. c£c , ) , 

(4.6) — ^ ^Af(0,£cEc t + t;(a,b)). 



Proof. The first assertion follows directly from Proposition 4.9 For the 
second assertion observe that by |Theorem 4. 13] the limiting distribution is the 
convolution of the normal distributions R(a, b) and the limiting distribution 
in ([45]). □ 

Corollary 4.15. In the situation of Theorem 4-lS\ assume that additionally 
iE(jW) -> b. 6 = Wi = 0, J^i = 0, ^ = | EI=i + 
E(^)-E(jf ) )) E(f), = Er=i E(J i (7V) )(E(jf ))+Er= 1 i ^ 
f = f(a, b) &e i/ie vector with coordinates fi = + + 6i— i + 26^ + 2 Yli=i a i- 
Then, 

t ( i ; 3/2 + e j ^ AT(0, If Ef + «(a, b)). 

Finally, we note for the convergence of the variance of R.W : 



Lemma 4.16. In the situation of Theorem J^.13 assume that additionally 
^E(jW) -> b ; and that ^J (Ar) is bounded. Then, 

Var^W) 

_ ^( a ,b). 

Proof. We have 

E((RW) 2 |jW) =Var(Jr tU (L( iV ) ) jW)) 

m 

= ^E(^+4S ) -4 JV) )^ + 42) 
i=i 
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and by our assumptions E((R.W) 2 | jW) /N 3 converges boundedly to v(a., b). 
Hence E((RW) 2 ) /iV 3 = E ( ( (RW ) 2 1 jW) ) /iV 3 — > u(a, b) . □ 



4.3. Unrestricted number of parts. We first consider the total number 
T(L)(g) of L-admissible partitions, without restrictions on the number of 



parts. In this case clearly A k -i+x(k) (as defined in Proposition 4.5) has 



a binomial distribution with parameters n = Lj. and p = j^rj, and hence 
each Jj can be represented as a sum of independent binomial variables. 
Furthermore, the covariance of A k _i + x(k) and Ak—j+x (fe) can be computed 
as 

Lk (.,.., ij 



Cov(A k _ i+1 (k),A k _ j+1 (k)) 
We therefore have 



fe + 1 



k + 1 



Lemma 4.17. Consider (Jj_, . . . , J m ) as defined in Proposition 4-5. Then, 
(4.7) 



E(J.) 

(4.8) Var(Ji) 

(4.9) Cov(Jj, Jj) 



. .fe+T' 

fc=i 



m T m j 



mm 



(fe + 1) 2 ' 



fe + 1 



E 



La 



, (fe + 1) 2 ' 

fe=max(i ,j) 



fc=max(j,j) 

Moreover, straightforward computations lead to the exact expectation 



value of the random variable defined in (4.1 ). We will need this for comparison 



to the basic specialization of Demazure modules in §4.5.2 Note that an 



Corollary 4.15 



asymptotic approximation (by Proposition 4.10) to this mean is used in 



Lemma 4.18. Consider the random variables defined in (4.1), and let 
^ = ET=i fcTT, U = Ek=\L k . Then, for E(T) = E(Y) + E(Q(L, J)) we 
have 



E(Y) 
E(Q(L,J)) 



1 



1 v — v I 



E-. 2 +E 



i=i 



4^ i + l 

i=l 



-Li 



m m m Y ' _|_ 2) 
£ ^ + £ iSiti + £ ^ 



i=l 



i=l 



-6(i + l)" 



Theorem 4.19. Let ^L^ -> a +- 0. Then, ^E(jf°) i^™ =i ^ /or 



eac/i i and 



JW - E(JW) _^ 

ATl/2 : 



AA(0,S) 
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where 



\j=mm(i,j) 

fc=max(i ,j) 



ak 
k + l 



fc=max(i ,j) 



ak 

+ 1) 2 ' 



Proof. The first assertion is clear. For the second assertion let B^ denote 
the underlying total occupancy statistic. In §4.1 it was shown that the 
components ~B(i)( N ^ are indendent multinomial distributions with parameters 
1 A! " wl >••' "<■' — — ■"' '"' — ^j- , and covariances L^E(i), 



o 



L- and u(z), where u(i] 

£(«) = diag(u(i)) - u*(*)u(i). By 
have 



u(i)i 



Theorem 4.7 



1 rW 



and since 



a,-, we 



B 



(JV) 



E(B™) 



#1/2 



(/V(0,aiS(l)),...,/V(0,a m E(m))) 



where the components on the right-hand side are independent. Since is 
a linear image of b[^ it is clear that J^' is asymptotically normal. The 
assertion about the covariance matrix is obvious. □ 



It is clear from Corollary 4.15 that under the conditions of Theorem 4.19 the 
number of admissible partitions is asymptotically normal, with expectation 
of order N 2 and variance of order jV 3 . The final condition of Proposition 4.10 
is fulfilled since the coordinates of the B(i) are binomially distributed and 
therefore have central fourth moments of order L 2 . The boundedness con- 



dition of 



Lemma 4.16 



explicit results for the case = (0, N) £ Z 



is fulfilled since < X^fcLi • We give some 



+ • 



(4.10) 
(4.11) 



E(TW) 
^VarCTW) 



, T ,m(4m + 5) . Nm(2m+1) 

iV ^ TV") + 



(4.12) ^ 



1 m 



12(m + 1) ' ' 12 m + 1 
m(m + 2) (2m + l)(2m + 3) 
45(m + l) 2 
m{m + 2) 
+ 36(m + l) 2 ' 

m(m + 2) 
6(m+ 1) ' 



Let us emphasize explicitly the representation-theoretic implications based 
on our discussion in[£ 



Theorem 4.20. Consider a sequence 



fC 2 )* L i *• • ■ *(C m+1 )* Lm ' of fusion modules of the current algebra sl2®C[t}. 



if Mm 



(JV) 



(JV) r (JV) 



a / 0, i/ien i/ie sequence of basic specializations 



Liffl),a)(q) behaves 



asymptotically normal with individual means 
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k=i k=l 

and variance, as N — >• oo ; 



k=i 1 

• oo ; 

1 



(4.14) ^^^ ) ^i fE f« + t,(a J b). 

Here, the vectors a, h, f, the function v, and the matrix £ are given as 
a = (ai,...,a m ) = lim — (L^\ . . . , 1^'), 

A— >oo iv 

m 

i-1 

/i = Of + 6f +1 + + 26^ + 2 a^, 

£=1 

j m 

u(a,b) = — } ](aj + b i+ i - &»)6i(aj + 6 i+ i) (Wiere 6 m+i = Oj, 



12 

i=l 



2j,- ,• = mm (i, i ) > ii > t ttt. 

M ^ fc + 1 J ^ (fc + 1) 2 

fc=max(ij) fe=max(i,j) 

Proof. This is simply a re- formulation of the results obtained in | Corollary 4.15 



and 
and 



Lemma 4.18[ For the convergence of the variance see Proposition 4.10 



Lemma 4.16[ □ 



4.4. T he central restricted case. We consider the "central region" (see 
§4.1. 2\ . Let sat = [E(S'^ Ar) )J = |_| E™1 iL i\- 11 is clear from the above that 



the underlying occupancy distribution is the conditional distribution 

m k 

B^ } = (Y(l), . . . , Y(m)) | £ Y ^ = s *> 

k=l i=0 

where Y(l), . . . , Y(m) are independent random vectors, each Y(/c) is multi- 
nomial with parameters L k N ^ and Po = ■ ■ ■ = Pk = jripi • This conditioning 
has the following effect on the asymptotic distribution. 

Theorem 4.21. Let m > 1, s N = [E(S^ } )\ and u = (u(l), u(2), u(mj) 
with u(k) = (jtj, . . . , j^rj)- Assume that -^L^) — > a^O and let c 2 (a) = 
h YJk=i K k + 2 ) a k- Then, 



'N 
where 



(4.15) Eij(k,£) = a k 5 k>e (^-^8 itj - ^^2 ) 
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(k-2i)(£-2j) 

~ aka U(k + i)(£ + iW(*y 

Proof. For ease of exposition we give the proof for the one component case 
LW = (l[ N \ . . . , L^} 1 , ) = (0, iV), the generalisation is straightfor- 
ward. Let := H(m)( N \ and X\, . . . , Xjy be i.i.d. random variables 
uniform on {0, . . . , m}. Let 5jy := 5Zj=i -^i> A 4 = m /2, <r 2 = j^m^m + 2). 
The probability generating function of B^-* is 

E(II & ) = l xSN l (*o + *!* + ...+ tmX m ) N /[(m + l) n P(S N = s N )j. 

i=0 

Hence the joint distribution is given by 



P(B, 



(JV) 



fed,...,S^0 =*;„») = L N k )/[(m + l) N P(S N = s N )) 



with the constraints that YliLo ki = N and X^iLi^i = s iV- Since there 
are two linearly independent linear constraints on the values of we 
expect a (m — l)-dimensional limiting distribution. Let xq, . . . ,x m be real 
numbers with YaLq x « = an d YhLq = 0> an d ^ ki = + y/Nxi. By 
Stirling's approximation for the factorials for the numerator and the local 
limit theorem for lattice distributions for the denominator we see 

(^""^(B^ = k) — > j==j V(m + l)m+i a 2 e - n ^ 1 (4+-+^), 
A check that the expression on the right-hand side is (considered as a function 



of X21 • • • , x m , say) the marginal density of (A/"(0, ^))2,...,m with S as ir|4.15 



concludes the proof. □ 

For the convergence of moments we have here 
Proposition 4.22. Under the conditions of \Theorem 4-21 



sty ->„„(*) am NH<tn^) 4l!M| 



Furthermore, E(B < (ib)W - E(B i (A:)( Ar ))) 4 /iV 3 — ► 0. 

Proof. Again we restrict the exposition to the one component case and use 
the same notation as in the proof of Theorem 4.21 From the generating 
function given there we get 

( iv) N P(S N -i = s N -i) 



(4.16) E(B 



m + 1 P(SV = sat) 



(4.17) E(5f)) 2 = E(i?f)) + iV(iV - 1)P( ^- 2 = SjV - 2z) 



(m + 1) 2 P(5jv = s jv ) 



(4.18) E(S W 5 W ^ - ^ " 1} P( - Sn - 2 = SN ~ i ~ ^ 



j 



(m + 1) 2 P(S N = s N ) 
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By the local central limit theorem for lattice distributions \17\ Corollary 
VIII. 3] we have 

1 (k-N^f / ,1/0, 

P(S N = k) = —e-^(l + 0(N-^ 

if | ^2JVo-^ I ^ Applying this to the numerator and denominator shows 

that the quotients q r (N) := P ^'p(g^ S j^ r )"' > are asymptotically of the form 

q r (N) = 1 - 2^]v + 0(N~ 3 / 2 ) wh ere c r = r 2 (m/2 - i) 2 . 

Now, the asympotic assertion about the expectation follows immediately 



from the local limit theorem, applied to numerator and denominator in 4.16 
The asymptotic assertion about the variance/covariance follows from the 
formulae above using the asymptotic form of qi(N),q2(N). Concerning the 
asserted convergence of the central fourth moment, note that the r — th 
factorial moment of B- is 

E(Bf)(Bf)-l)...(Bf)-r + l)) 

N(N - 1) ... (N - r + 1) P(Sjv_ r = s N - ri) 



(m + l) r P(S N = s N ) 

After expressing the central fourth moment as a linear combination of factorial 
moments and plugging in the asymptotical expressions for the q r (N), a little 
algebra yields 

E(£f° - E(flf°) 4 = iV 3 (c 4 + 6c 2 - 4c 3 - 4ci) + 0(iV 5/2 ) = 0(iV 5 / 2 ) 
since C4 + 6C2 — 4c3 — 4ci = 0. □ 
Remark 4.23. Let — ► a ^ 0. A comparison to the unrestricted 



case, discussed in |§4.3[ shows that asymptotically the underlying total 
occupancy distributions are quite similar. They concentrate around the 
same expectations. In the unrestricted case the components of the limiting 

distribution -^=(B^ — u) — Z are independent normal vectors with 

Cov(Z(fc)) = a fc (diag(u(A;)) - u(fc)*u(Jfe). 

The components stay normal in the central restricted case, but the restriction 
causes an additional negative correlation 

n I7(u\ view a k a e (2k - i)(2£ - j) 
Cov(Z(k)i,Z(£)j) - 



4(fc + l)(Z + l)c7 2 (a) 



between the components. This in turn forces the elements of the asymptotic 
covariance £ of to be smaller than in the unrestricted case, we compute 

j, restricted — unrestricted 2a'-*(a) < "(^)'"(^) 

with c(i) = YZ=i 
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Since is a linear ima ge of its d istri bution is also as ymptotically 

""that is 



normal and it is clear from 



Theorem 4.21 



and 



Corollary 4.15 



asymptotically normal and the preceding results show that the expectation 
resp. variance are of N resp. N 3 , but the variance in the restricted 

case will (on the TV 3 scale) be smaller than in the unrestricted case. 

Again, let us emphasize the implications for fusion modules of the current 
algebra sfe ® C[t]. 

Theorem 4.24. Consider a sequence J-(b^ , L^' . 



,L^°) = (C 1 )"'i * 

(C 2 )* L i *• • •*(C m+1 )* im > of fusion modules of the current algebra sh®G\t}. 
If ^(L±\ . . . , Lm ^) —J- a 7^ 0, then the sequence of central string functions 

T((l[ N \ lJ£P), sn )(q), i-e. sjy = [\YliLi^i\ ' behaves asymptotically 
normal with asymptotic mean 

(4.19) 



■il\*b 



(AT) 



E 

i=l 



m 

E Qfc 
fe+1 

k=i 



+ 1 



E 

k=i 



ak 
fe + 1 



i-1 



Yj ak 



k=l 



and asymptotic variance 



(4.20) 



1 J 



ifEf* + w(a,b). 



Here, the vectors a, b, f , the function v, and the matrix £ are given as 
a= (ai,...,a m ) = lim -^-(l[ N \ . . . , Lffl), 

TV— >oo iv 



E 

£=i 



ai 



i-1 



fi 

v(&, b) 



H + b i+ i + bi-i + 2bi + 2 ^2 at, 

i=i 

. m 

— y_ '(oj + 6 i+ i - bi)bi(ai + (Wiere 6. 



12 



Yn+l 



i=i 



mm («,j) ^2 

fc=max(i j) 



«-A: 



2<7 2 ( 



fc + 1 



'3 



E 



fc=max(i ,j) 



(fe + 1) 



where a 2 (a) 



1 ZX=i + 2 W> and c(i) 



12 



y^m fc+l-i 



Theorem 4.21 



Proof. This is again a re-formulation of 
ear image of B^), the asserted convergences follow from Proposition 4.10 



Since is a lin- 



Lemma 4.16, Proposition 4.22, and Remark 4.23, Note that by Proposi- 



tion 4.10 the leading term of the expectation values, i.e. the coefficient of 
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N 2 , depend only on a and b in the same manner as in the unrestricted case. 
Therefore, to derive the asymptotic mean given in (4.19) one simply has to 
replace the in the quadratic terms of (4.13) by their limit values a^. □ 



(l[ N \ . . . , i^2 1; Lm^) = (0, N) to the corresponding unrestricted one. We 



Let us explicitly compare the central restricted one component case L/^ 

(TV) 

'1 , • ■ 

rewrite 

m 

rW\2 



Q( L W JW) = £(JJ 

i=l 

m 

= £(jf> - E( J^)) 2 + 2 - E(JD) E(jf> ) 

=i 

m, 

+E( E (^)) 2 ' 



i=l i=l 

m, 

i=l 

and let in the sequel 



G W = ^(jf)-E(jf))) 2 , 

i=l 

m 

AfW = 2 jf } - E(jf°)) E( J f (7V ; 



It is straightforward to compute that in the case of jointly normal variables 
jf° we have Var( E™i(#°) 2 ) = Var(GW) + Var(MW). 



Lemma 4.25. Lei B = N(Nu, NT,) with £ given as in Theorem ^.7, Then, 
G^ has expectation 

EC« = V Vta<jf>> = jV(m-l)(m» + Mro + 6) 
^ v ? y 15 (m + 1) m + 2) 

and variance 

(7V) _ iV 2 (m - l)(llm 5 + 99m 4 + 16482m 3 - 2372m 2 + 6696m + 2304) 
' ~ 6300(m + l) 2 (m + 2) 2 ' 

Furthermore, the variance of is given by 

(Ar) _ N 3 m(m - l)(m 3 + 247m 2 + 136m + 12) 
' ~ 180 (m+ l) 2 (m + 2) ~' 

Thus M^ N > is of order of magnitude iV 3//2 m while G^ N ' is of order of 
magnitude mN . Therefore the positive random variable G^ will cause a no- 
ticeable distortion of the asymptotic normal distribution of [Q(L^ N \j^ N ^) — 
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E(Q(LW, J w )))/JY 3 / 2 unless iV » y/m. Consequently, one has 
(4.21) 
ET (7V) 



9 m(4m + 5) AT 



in 



m 3 + 19m 2 — 14m — 6 A^ m — 1 



12 m + 1 



15 m + 1 



m + 2 



4 m + 1 



Remark 4.26. (1) We computed the expectation of T^ N ' and the variance 
of Q(Ij( n \ j( N ^) here under the simplifying assumption that the occupancy 
distribution is normal. Note that we know from Proposition 4.10 and 



Theorem 4.21 that the leading terms in N are correct. 



(2) Comparing these results to (4.10) given at the end of the previous 



section, we see that the additional restriction has only small effects on the 
asymptotic distribution of . Expectation and shape remain essentially 
unaltered, only the variance has decreased. This was to be anticipated from 
the probabilistic setup: the conditioning is "in line" with the law of large 
numbers. 

(3) The cases j^S^n) — > e' 7^ e:= \ Y^k=\ can be treated by standard 
large deviation techniques. Again we may safely expect to find normality of 
the asymptotic distributions. 

(4) In terms of fusion modules the results in Lemma 4.25 and ( |4.21 1 apply 
to the central string function f((0,AQ [Nm/2\) in T jO^N) = (C m )* N and 
can be re- formulated accordingly (see Theorem 4.24[ ). 



4.5. Comparison to previously known results. We compare our find- 
ings to established results in the literature 0E1E1E2]. 



4.5.1. The unrestricted one component case. Consider first the one component 
case L = (l[ N \ . . . , L^J_^_, Lm ) = (0, N). In this case the distribution of 
has been investigated under the heading "generalized Galois numbers" 
by Bliem and Kousidis [5] and later by Janson [22]. These authors studied 
the random variables with probability generating function 



(4.22) 



E(g 



1 V 

( m + l)tf 1^ 

(fc ,...,fc m )eN™ +1 

k +...+k m =N 



Note that if we let Bq = L m —J m , B\ = J m — J m -i, 



N 

ko, . . . , k m 
■ ■ , B m -i = 



Ji—J\ , B Tl 



J\ it becomes evident that this distribution coincides with the distribution 
of Y( N ) in the one component case. Now, Bliem and Kousidis showed 

Theorem 4.27 ((5J Theorem 3.5]). Consider the random variables 7^ 
defined through (4.22). Then, 



E(yW) 

Var(yW) 



1 



m 



N(N-l), 



4 m + 1 
1 (m+1) 2 
72 (m+1) 2 



1 



N(N-l)(2N + 5), 
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and 



y( n ) - E(y( Ar )) 



M(0, 



1 (m + 1) 2 
36 (m + 1) 



Janson derived the same result in a variety of ways, proved a corresponding 
local central limit theorem, and gave an overview of different interpretations 
of the distribution of . In addition he showed joint convergence of yw 
and B^: 

Theorem 4.28 ([221 Theorem 2.4]). Let B^, yW 6e as above. Then, 

/yW-E(yW) bW-e(bW) 



(4.23) 



AT3/2 
d 



AT(0, 



JVV2 
1 (m + 1) 2 



1 



),AT(0,E) , 



36 (m+1) 2 

where the constituents on the right hand side are independent, and the matrix 
£ is given by = (^i,j - — ^ * 



m+1 - 



Let us compare these results to our findings above. We obtain from 
ILemma 4.181 

e(yW) = 1 -^-n*- 1 -^-n, 



4m + 1 



4m + 1 



which agrees with the expectation given in Theorem 4.27 Further, we have 
N i 



and find that 



m+1' 



a 



m + 1 

m(m + 1) 



a m -i — 0, a r 



1, 



iV 2 



and 



w(a, b) 



12(m 



1 m(m + 2) 
36 (m + 1) 2 



Finally, c = 0, where c = c(a, b) is as in Corollary 4.14. Thus, by Corol- 
lary 4.14 we have 

m(m + 2) 
JV3/2 '36 (m + 1) 2 j ' 

which is equivalent to the weak convergence assertion in |Theorem 4.27 and 



hence establishes an independent proof. Moreover, by Corollary 4.14 
e (LW,j(*0)-d(*> d 



iV 3 / 2 



that together with Theorem 4.1 independently proves Janson's Theorem 4.28 



Note that Theorem 4.28 as it stands does not generalize to more general 
distributions. As an example let be multinomial with parameters N, p 
where p is not uniform. Here we get from Corollary 4.14 that 

e(LW.jm-dW A 
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where i>i(p) = \ { YliLoPi ~ (YliLoPi ) 2 ) • T he corresponding joint limiting 
distribution (on the right hand side of (4.23)) is normal, but the constituents 
are not independent. 



4.5.2. The unrestricted two component case. For the two component case of 



L = (L 1 , 
we find 

E(T L ) 



For K 



i.e. L m = M, Lfc = K for a k < m, and all other Lj = 0, 



1 m(4m + 5) ~ 1 m(2m + 1) , , 1 kUk + 5) o 
' M H -M+— — K 



12 m + 1 

1 fc(2£; + l 



12 m + 1 



12 jfe+1 



+: 



E(T L ) 



12 fc + 1 
1 this simplifies to 

1 m(4m + 5) 



if + l mi f M + - fc(fc + 2) KM. 
2 6 m + 1 



o 1 m(8m + 7) + 2k(k + 2) „ r jb 
M H M + 



12 m + 1 12 m + 1 2 

Let us compare this to the only known result for a two component case in 
the literature, the case of Demazure modules V^(A) associated to the affine 
Kac-Moody algebra s^. Fix the highest weight A = mAo+nAi. The random 
variables X w having probability generating function the basic specialization 



of the character x(K«(A)) are given due to Proposition 3.1 by translations 
and rotations (averaging over the random variable Sl^) a s follows 



(4.24) 
(4.25) 



Equations (3.4) and (3.5) read as: 

X w = N 2 m + N(N + l)n + (-2N - 1) • S Lw + T u 
Equations (3.6) and (3.7) read as: 

l)n 



N 2 m + N(N 



2N-S Lw +T u 



The four cases covered here correspond to the four cases in Theorem 4.1]. 
Let us restrict for simplicity reasons to (4.25) for w = (sxSq) , and compare 



our findings to [TJ Theorem 4.1], where the corresponding case is (4.1)] 
for even N. 

Theorem 4.29 (j7l (4.1) in Theorem 4.1]). For L with entries except 
L m = 1, L m + n = 2N — 1, and with U = 2N — 1, u = m + n one has 



(4.26) 



E(N 2 m + N(N - l)n - 2S L N + T L ) 



2Um(m + 2) + U(U 



12(u + 1) 



l)u(u + 2) U 



— 1 u m 
2~2 + ~2' 



We can re-establish (4.26) by the computation of the left-hand side through 



the linearity of E(.) and the mean of the random variables Sl, 7l. 

Proof. From the above we have in this case 

iu + 7) + 2m(m + 2 



J_n(4, + 5) 
v ' 12 u + 1 



1 u( 
12~ 



u + l 



U + 



in 
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and E(S'l) = \ {mL m + uL u ) = \{m + uU) = ^(2uN — n), which gives 

(J J + "H 2 

E(iV 2 m + N(N - l)n - 2S L iV) = -N 2 u = - ' u. 

It is now readily seen that both results agree. □ 

Since the results in [7] are derived differently and independently from our 
arguments through the analysis of Demazure's recursive character formula, 
this nicely confirms both, the correctness of the theory given here and of 
Theorem 4.1]. 

Finally, we settle a question that was posed by Bliem and Kousidis in [6]. 

Lemma 4.30 (Cf. [6] Conjecture 8.3]). Fix a dominant integral weight A 
and a sequence (w^) in W &s such that l(w^) — > oo. Let be the joint 
distribution of the degree and the finite weight in V w (k) (A) . Let fi^ be the 
distribution obtained from by normalizing to a probability distribution 
and rescaling the two coordinates individually so that supp(fl^) just fits into 
the rectangle [0, 1] x [—1, 1]. Then, as k — > oo, 

~(k) w x 

V > Of <e ,A)+2 \, 

where c = Qq + a\ denotes the canonical central element. 

Proof. We consider only the Demazure modules V( SlSo );v(mAo + nAi) as the 
other cases can be derived equivalently. Let d^ denote the maximal degree in 
these Demazure modules, i.e. d^ = N 2 m+N(N —l)n. Let u = m+n = (c, A) 
denote the level of the representation, and consider the random variable 
with probability generating function given by the basic specialization of our 
Demazure module, that is 

X N = d N - 2NS UN) + TW 

= B(d N - 2NE(S UN) ) + rW) - 2N(S UN) - E(5 L(W) )) 

+ ( T w_ E (tW)). 

The probability distribution of Xjv and j^Xjy is the first coordinate of fj,^ N ' 

and pS N } for the Weyl group element = (siso) N , respectively. Now, 

equivalent to the asserted weak convergence of jl^ N ^ we have 

X N d u+2 

d^ 3(it + 1) ' 

since 

(N) + T( 7V )) 1 4« + 5 u + 2 

ahi > 3 u + 1 ~ 3(u + 1)' 



and by (4.3) and Corollary 4.15 we have the convergences in distribution 
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Since it is already known that the second coordinate of concentrates in 
by the factorization phenomenon (3.10), the claim follows. □ 



The basic specialization can be interpreted as the generating function of 
dimensions of si2 representations i n V W (A) splitt ed by the degree, i.e. the 



coefficient k in e ka °. Therefore, Lemma 4.30 depicts their asymptotic 



concentration along the degree as a rational expression in the level of the 
Demazure module. 

One can easily compute that the asymptotic concentration in the central 
restricted two component case (i.e. along a central string function) is the 
same. This exhibits the dominating weights in Demazure modules. Note that 
this is not deducible from [61 [7] due to the limitations of their arguments 
imposed by the non-positivity of Demazure's character formula. The basic 
argument that the concentrations coincide is that the expectation values of 
the random variables in the central restricted two component case have the 
same leading terms (in N) as the unrestricted ones. 

Let us briefly describe the central string functions in fusion modules whose 
study would generalize this asymptotic concentration. Let A^ 1 £ z™ xm 
be the matrix with A~j = min(i,j), and consider the vector £ = j4 _1 L. In 



particular, t\ 
and following 



^2 Li and t m = ^ihi. Then by (4.2) we have E(S*l 



14.1.2 



1/ 

2 'mj 



the g-supernomials T(L, \lm) describe the central string 
functions in fusion modules (resp. T(L, \l m =L |)), that carry the dominant 
weight multiplicities. Due to |35| Definition 2.3] and |35[ (3.1)] the translation 



-,L t A~ 1 L 



(4.27) 



n 1 1 If: 



T(L, ^^r 



i=0 



is a polynomial with degree at most \L l A 1 L. The reader is referred to 
Figure 2 for an example in the case L = (0, 4), where we have q~ 8 T((0, 4), 4) = 



l + q + 3q 2 + 3q 3 + 4q 4 + 3g 5 + 2q 6 + q 7 + q 8 . 



5. Local central limit theorems 



It should be possible to prove the following local central limit theorem 
along the lines of [22\ §6]. We pose them here as conjectures. 



Conjecture 5.1. In the notation of Theorem 4-20 let Xp N = X (n) 

denote a random variable with probability generating function the normalized 
basic specialization of the fusion module J-(b^\ h\ , . . . , Lm ^). Denote its 



mean /^n - 
ink as N 

(5.1) 



-> oo, 



(iv) and variance air 



a 



(AT) 



r (N) 



Then, uniformly 



^a N ■ P(X Tn = k) = e 



-{k- l i N )/2<j 2 N 



+ 0(1). 



:!2 
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Here, a'jq can be replaced by the explicit expression A^ 3 (|fSf' + t>(a, b)) /? 



rom 



(4.20). In particular, the dimension of the sfo submodule in J~n of degree k 
grows as (5.1 ). 

The analogue for the central restricted case is: 



Conjecture 5.2. In the notation of Theorem 4-24 and (4.27) let Sn denote 
a random variable with probability generating function the normalized central 
string function 



1,(7V),(A0 , , 

- 2 'l l m TClS > 



,Sn)- 



Assume ^ a ^ 0, then ±l[ N) 
we let 



1_,(A) 
N<"in 



/' 

2 



lim 1 



.,L^ m '),s N 



2^1 



lim 



A 3 



(7. 



(JV) 



A— >oo 

Then, uniformly in k as N — > oo, 
(5.2) V^ko- ■ P(S N = k) 



-(fc-/i)/2cr 2 



+ o(l) 



In particular, the dimension of the t) C sfe weight space with coordinates 
\lm^cti and —k5 grows as (5.2). 

Note that those statements are the first hints towards an asymptotic 
description of the dimensions of fundamental submodules in fusion products 
available in the literature. 

The proofs should work verbatim as the proof of |22l Theorem 1.5] once 
one establishes an estimate for the generating function T(L)(q)/T(L)(l) ana- 
loguous to [22l Lemma 6.1], and for the g-supernomials T(L, a)(q)/T(L, 
analoguous to [22], (6. 2), (6. 3)], respectively. 



6. Fusion of symmetric power representations 
The Kostka numbers are the coefficients in the expansion 

(6.i) nM x )=E^-^ x ) 

of the product of complete symmetric functions h^ in terms of the Schur 
functions s v - The Kostka polynomials K Vyfl (q) generalize the Kostka numbers 
in the sense that K V)fl (l) = K Vyfl . They give the transition matrix between 
the Schur function and Hall-Littlewood function P^, i.e. 

(6-2) ^(x)=^^(g).P M (x,g). 
A standard reference for the above functions is [29) . 
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Now, the g-supernomial 5^ i M (g) [21] |25| 134] 1351 136] is defined as the 
combination of (6.1) and (6.2), i.e. as the transition between the above 
product of complete symmetric functions and Hall-Littlewood functions 

(6-3) SeM = Yl Kr >£ " K vAl)- 

An explicit form of S^^(q) is proven in [21], Proposition 5.1], where fi = 
Qui, . . . , fjL m ) is a partition and £ G Z™ a composition such that \fj,\ = |£| = M, 
as 



3.4) 



with 



^)=e^ m) n 



l<o<n-l 
1< i</J,i 



(a). 

(a) 
4-1 • 



n— 1 /ii 



(a+l) (a) 



a=0 i=l 

and where the sum ^ s indexed over the sequences of Young diagrams 

i/W . . , i/ 71 - 1 ) such that 



C i/W C i/« C • • • C zA- 1 ) 
»| = & + ... + £, for 1 <a<n-l. 



Remark 6.1. For n = 2 and arbitrary yu, (6.3) agrees with the definition of 
g-supernomials as given by Schilling and Warnaar [35]. See |25[ §3.1] for a 
detailed discussion. 

We define a slight variant of the above g-supernomials, which describes 
the string functions in the fusion product of sl r +i symmetric power represen- 
tations F„, = V r MlWl * V^ 2U11 * ... * Vfj, mU11 . That is, 



where for the partition /j = Qui, . . . , (J, m ) we set n(/i) (Cf. [Ml (3.10)], [26 
§2.1], |35[ §2.1]) to be the normalization constant 



3.6) 



n 



Qu) = - = min (Mi> Mi)- 

l<i<j'<m 



i=l 



Note that this normalization ensures that q n (v) K„ a(q 1 ) is a polynomial in 



Then, we have a fermionic formula (a positive sum of products of q- 
binomial coefficients) for the graded character of the above fusion product 



Ffj,. That is, 
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Proposition 6.2. Let \i = (m, . . . ,fi m ) be a partition of M. Then, 

£ weight 

Proof. Let denote the monomial symmetric functions. Then, with 

K V M = Q^KvAl' 1 ) where n(n) = £ f (i - 1)W a s m M (3-10)] one 
has [24, Corollary 7.6]: 

X{^n) = Y ' k*iAq) 
= s v - k v Aq) 

T] \-M 

= J2 { Y K v,r m q -KriM 

7] \-M ^ £ weight ' 

= Y ( Y K v,rKr,Ai)j - mi 

f weight ^ r] \-M ' 

= Y • m « 

f weight 

Note that all partitions except /j, have at most r entries, corresponding to 
the rank of the Lie algebra. □ 

Remark 6.3. For the graded character of fusion of fundamental representa- 
tions *jV(uii ), Chari and Loktev prove an equivalent fermionic formula [9, 
Proposition 2.1.4]. 

Kirillov |25| [26] is a great source of various combinatorial, geometric and 
statistical interpretations of q-supernomials S^u(q). Let us shortly remark 
on the geometric one. 

Remark 6.4 (Cf. [3 [^[271 SB]). As pointed out by Kirillov [25J §1.4] it has 
been proven by Shimomura |37] that the g-supernomials count the number 
of rational points Fl^ (F 9 ) over the finite field F q of the unipotent partial 
flag variety Fl^. To be precise, for a composition £ G of re, a £-flag 
in a n-dimensional vector space V is a sequence V% C • • • C V r such that 
dim Vi = £l + • • • + The set of all such flags is the partial flag variety Fl^. 
We let Fl^ C Fl^ be the subset of the partial flag variety Fl^ consisting of the 
set of all £-flags F G FL- fixed by a unipotent endomorphism u G G1(V) of 
type fj, (a partition of n that describes the Jordan canonical form of u). Then, 
Fljf is a closed subvariety of Fl^ , the so-called unipotent partial flag variety. 
Now, Shimomura [37] proves that the q-supernomials count the number of 



Fg-rational points in F!|. That is, with n(/i) as in (6.6) one has 
(6.7) #Fr|(F,) = fMStoiq- 1 ) = S^q). 
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In particular, the basic specialization of the fusion module gives the 
number of Fg-rational points in JJ^ Fl^ , 



(6.8) 



x(J>)(?) 



E 



Our Proposition 6.2 exhibits the objects that have to be analyzed in 
order to establish a general central limit theorem along the same lines as 



Theorem 4.20 The explicit expression (6.5) shows that one can interpret the 



g-supernomials again as mixtures of probability distributions. For an Ansatz 
let 



K, 



Here, f^nil) would immitate the inversion statistic, and X„ £ the mixture 
distribution. It should be straightforward to check the reductions to the 
distributions investigated in §4 in the case of g-supernomials as defined by 
Schilling and Warnaar (see Remark 6.1). We pose a conjecture for further 
research. 

Conjecture 6.5. Consider the sequence of type A fusion modules of sym- 
metric power representations 



(AT) 



* V 2 * * v 1 



associated to the partition 
i 's. Assume that as N — > oo we have 



(i/r,2/ 2 (JV \ 



lh ) with f- N ^ -many 



1 

N 



(/i 



(JV) AN) 
5 J 2 ' 



a / 0. 



Then, the central string functions and the basic specialization of T„(n) behaves 
asymptotically normal as N — > oo. 
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